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If cosmic magnetic fields are indeed produced during inflation, they are likely to be correlated 
with the scalar metric perturbations that are responsible for the Cosmic Microwave Background 
anisotropics and Large Scale Structure. Within an archetypical model of inflationary magneto- 
genesis, we show that there exists a new simple consistency relation for the non-Gaussian cross 
correlation function of the scalar metric perturbation with two powers of the magnetic field in the 
squeezed limit where the momentum of the metric perturbation vanishes. We emphasize that such 
a consistency relation turns out to be extremely useful to test some recent calculations in the lit- 
erature. Apart from primordial non-Gaussianity induced by the curvature perturbations, such a 
cross correlation might provide a new observational probe of inflation and can in principle reveal 
the primordial nature of cosmic magnetic fields. 
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Cosmological inflation has become a very successful 
and dominant paradigm to understand the initial con- 
ditions for the Cosmic Microwave Background (CMB) 
anisotropies and Large Scale Structure formation [1]. In 
particular, the precise observations by the Wilkinson Mi- 
crowave Anisotropy Probe have provided a tremendous 
support for the inflationary paradigm 0]. These obser- 
vations are remarkably consistent with the predictions of 
minimal single held slow-roll inflationary models and in 
order to test further, it has become really important to 
look for new observables, such as non-Gaussianity that 
can be used to discriminate between different models of 
inflation and possibly rule out the minimal model. The 
detection of non-Gaussianity has therefore become one of 
the primary goals of present and future experiments Q . 
But it is very model dependent, whether primordial non- 
Gaussianity from inflation will be observable. 

Another possibility which has been contemplated is 
that the inflaton, or some other light field during infla- 
tion, for example, the curvaton |4| is coupled directly to 
electromagnctism. Such a coupling breaks the conformal 
invariance of electromagnctism and as a result, large scale 
magnetic fields can be produced during inflation ||. It 
has been speculated that this could even provide the seed 
fields required for the galactic dynamo [6j and also lead 
to acoustic signatures in the CMB 0, Hj]. However, there 
exists a theoretical problem for such models to succeed, 
called the strong coupling problem and as discussed 
in [l(| , it is very difficult to avoid it while retaining gauge 
invariance unless one gives up on the gauge field being 
the actual electromagnetic field, and instead thinks of 
it as some novel hidden sector isocurvature gauge field. 
Recently, it has also been speculated that the strong cou- 
pling problem can be circumvented, if gauge invariance 
is broken in the UV in the effective four dimensional La- 
grangian, but manifest in five dimensions (llj . One might 



also worry that the usual "77-problem" of inflation could 
be aggravated by such a coupling of the inflaton to the 
gauge held [1 21 ] . Independent of whether any specific pro- 
posal to overcome the strong coupling problem may or 
may not work, the model is an archetypical model of mag- 
netogenesis, and as we discuss below, our motivations to 
study such a scenario are more general. 

In a model with a direct coupling between electro- 
magnetism and inflaton, the metric fluctuations will be 
correlated with the produced large scale magnetic fields. 
Apart from primordial non-Gaussianity, the three-point 
correlation of the metric perturbation with two powers 
of the magnetic field might provide a new observational 
probe of inflation that could also shed light on the pri- 
mordial nature of cosmic magnetic fields and has been 
studied recently in 
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13]. In this letter, we show that 
in the squeezed limit where the momentum of the scalar 
metric perturbation vanishes, there exists a simple con- 
sistency relation that can be used to obtain the three- 
point correlation function of the metric perturbation with 
two powers of the magnetic field in a simple way which 
can be used to check the results in the literature. 

In order to evaluate the cross correlation of the cur- 
vature perturbation with the magnetic fields, it is con- 
venient to define the magnetic non-linearity parameter 
&jvXj in terms of the correlation function of the curvature 
perturbation with the magnetic fields as 

(C(ri ) k 1 )B(rr,k 2 )-B(rr,k 3 )) 

= 6 AfL (2^) 3 <5( 3 )(k 1 + k 2 + k 3 )F c (fci)P B (fc 2 ) (1) 

where Pq and Pb are the power spectra of the comoving 
curvature perturbation and the magnetic fields, respec- 
tively and are defined as 

<C(r,k)C(r,k')) - (27r) 3 ^(k + k')P<(fc), (2) 
(B(r,k) -B(r,k')) = (2^) 3 ^ 3 )(k + k')P B (k). (3) 
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The conformal time r is defined by adr = dt where a(t) 
is the scale factor of the Friedmann-Lemaitre- Robertson- 
Walker metric ds 2 = —dt 2 + a 2 (t)dx 2 and tj denotes the 
conformal time at the end of inflation. 

The time-dependent coupling of the electromagnetic 
field to the background, can be parametrized by a cou- 
pling of the form X((f)F llv F tlv , where = d^A v -d v A ll 
is the electromagnetic field strength and the time depen- 
dence of the coupling is parametrized by its dependence 
on a slowly rolling background scalar field <f>, which we 
think of as being the inflaton for simplicity. When 
is momentum independent, it corresponds to a "local" 
non-linearity which can be obtained from the relation 



B = B< G > + -^C (G) B( G > 



(4) 



where B^ G ) and are the Gaussian fields. One can 
estimate the size of b^L by noting that the interac- 
tion Lagrangian between the scalar field and the elec- 
tromagnetic field is Cqbb oc \(4>)F 2 . By Taylor ex- 
panding the coupling in the inflaton fluctuations, \(<fi) = 
\{4> c ) + d<p\(4> c )54>, one obtains that the linear coupling 
between the inflaton fluctuation and the electromagnetic 
field Cqbb oc d ( j ) X{(f))5(j)F 2 . It is useful to express the 
scalar perturbations in terms of the comoving curvature 
perturbation £ which can be considered as the scalar per- 
turbation of the metric 



ds 2 = -dt 2 + a 2 (t) e 2 « 4 ' x W 



(5) 



on large scales where the time derivate of C vanishes. 
The comoving curvature perturbation is related to the 
inflation fluctuation by 5(f) — V2e( where the slow-roll 
parameter e is given by v2e = —<j>/H. With these defi- 
nitions, we have 



(6) 



which leads to Cqbb oc —\((/H)F 2 . In analogy with 
the analysis of [14| , we can compare it with the quadratic 

1/2 

teim £bb- The ratio is CqbbI^bb oc X/(HX) -P^ , and 
we would therefore expect 



b NL ~ 0(X/(HX)) 



(7) 



We will show later that in the squeezed limit, where 
the momentum of the curvature perturbation vanishes 
i.e. fci <C ^2,^3 = k, the magnetic consistency relation 
in fact gives 



? local 

"nl 



1_X 
H X 



(8) 



which agrees with the naive expectations. 

Despite a basic difference (which will be evident later), 
due to the conformal nature of electromagnetism when 



the coupling A is trivial, the consistency relation shares 
some familiarity with other consistency relations, which 
have previously proven to be extremely powerful in the 
literature. In flil . Il6l ] , a consistency relation of the bis- 
pectrum of the curvature perturbation was obtained, and 
in [14| , a new consistency check of the tree-level exchange 
diagrams, where a scalar or a graviton is propagating be- 
tween two pairs of external legs, was derived. Some of 
these consistency relations have been extended to higher 
orders using arguments related to conformal symmetries 
in [It], EH ■ Finally, the semiclassical relation for loop di- 
agrams in [Tj| |Ti(| made it possible to clarify the nature 
of IR loop effects in de Sitter. Some of these approaches 
are reviewed in the appendix of 19]. 

In the Coulomb gauge with Aq = and diA % = 0, the 
quadratic action for the electromagnetic vector field Ai 
becomes 



s,. 



d 4 x^X(<l>)F^F^ 



= \J d 3 xdrX^) (a? - ±(0^ djAif^j (9) 

Since the magnetic field is a divergence free vector field, 
the two-point correlation function of the magnetic field 
can be written as 

( Bl (r^)B\rX)) = ^ ^ - (Mt^MtX)) 

(10) 

and the magnetic field power spectrum becomes 

P B (fc) = ^(A(r,k).A(r,-k)). (II) 

By parametrizing the time-dependence of the coupling 
function as 



A(</>(r)) = Xi(t/t!) 



-2n 



(12) 



one finds that the spectral index, rig, of the magnetic 
field power spectrum is given by tib = 4 — 2n for n > 0. 

We are therefore interested in computing a basic corre- 
lation function such as (C(tij ki)Aj(Tj, \s.2)Aj[ri , ks)) in 
the squeezed limit k\ <C &2, &3- Since £ is frozen outside 
the horizon, the effect of in the squeezed limit is to 
locally rescale the background when computing the cor- 
relation functions on shorter scales given by ki , k% and 
one can write [TJI, [Tj| 

lim (C(Tj,k.i)Ai(Ti, k 2 )Aj(Ti, k 3 ) 

ki— >0 

= (C(r /) k 1 )(A i (T J ,k 2 )^(r J ,k 3 ) B >, (13) 

where (Ai(ri, k2)A,(T/, ks) B is the correlation function 
of the short wavelength modes due to the variation in the 
background produced by the long wavelength mode of (. 

In the absence of the coupling function A, the confor- 
mal invariance of the gauge field implies that it only feels 
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the background expansion through A which subsequently 
depends on the scale factor. Since Ai(r, k) doesn't feel 
the background for a trivial A and the only effect of 
£(t, k) in the squeezed limit is to locally rescale the back- 
ground as a — > as = e C ' B a, we expect that the correlation 
function (C(tz~> \ci)Ai(ri, k 2 )Aj(T/, k 3 )) will vanish in this 
limit to the leading order for a trivial A. 

To compute the correlation function for a non-trivial 
A, we need to write Ai(r, k) in terms of the Gaussian 
one with a trivial A. If we expand A^ = A(as) around a 
homogenous background value Aq = A(a), we get 



As = A + S\na- 
ama 



\ dX ° r 

Ao + -r, Cb 

a ma 



(14) 



Notice the resemblance with a 5N expansion. Using N = 
In a, we could also have written the above as A = Ao + 
X' N 6N + . . . , where \' N = dX/dN. 

Defining the linear Gaussian part of Ai to be A\ G ' and 
in the Coulomb gauge with a pump field S 2 — Ao, we 
can define a linear Gaussian canonical vector potential 
Uj = S(t)A[ G \ such that the quadratic action in © 
takes the form similar to a canonical scalar field with 
an effective time dependent mass term S" / S and can be 
written as 



S v 



drdrx 



S" 



(15) 



In the squeezed limit, all the effect of (b on the gauge 
field is captured by (fl"4"|) and therefore, we have 



(A l (T,x 2 )A j (T,x 3 )) B = ( — « i (r,X2)«j(T,x 3 ) 



— (vi(T, x 2 )uj(r, x 3 )) 

^0 



1 dX 
Ag d In a 



Cb (^(t,x 2 )uj(t,x 3 )) (16) 



which leads to 



(A{(t, x 2 )A # (t, x 3 )) b 
1 dX 



Aq d In a 



~ (A 4 (t,x 2 )A j (t,x 3 )) 
CB^i(r,x 2 )^(r,x 3 )) (17) 



i dx r d 3 k B 



and a Fourier transformation gives 

{A l {rM)A 3 {rM)) a 
-C(r,k s )(A 2 (r,k 2 )^(r,k 3 )) (18) 

Using dX/d\na = X/H. we then find from (fT3|) the 
squeezed limit consistency relation for the gauge field 

lim (C(r/, k^Aiin, k 2 )^(rj, k 3 )) 

fei — >0 



1 A 
HX 



(C(T I ,k 1 )C(Tl,-kl)) (MT I ,k 2 )A J (T I ,k 3 )) 

(19) 



Finally, by using the relation (flO|) , we find the consistency 
relation for the magnetic field to be 

(C(7j,ki)B(7j,ka).B(rr,k3)) 



k 2 + k 3 )P c (A; 1 )P B (fc 2 ) (20) 



= -l^(2^ 3 )(k 1 

in agreement with © . With the parametrization in (fTz 
we obtain the consistency relation 6jvl = — 4. 

In the squeezed limit, the consistency relation is quite 
general as an explicit form of the coupling function has 
not been used. But as we argue below, the approxima- 
tion used to obtain this consistency relation might only 
be trusted for n > 1 in (|12l) . To see this, note that for 
a canonical massless scalar field in de Sitter space, the 
pump field in (|15|) can be identified with the scale fac- 
tor, and one would have S"/S — a" /a with a oc 1/r. In 
de Sitter, there exists an apparent future event horizon 
and therefore, the long and the short wavelength modes 
decouple on super horizon scales. As a result, all the ef- 
fects of the long wavelength modes can be captured by 
their rescaling of the background of the short wavelength 
modes. In the more general case above, if we interpret 
S as an effective scale factor for the background experi- 
enced by Vi, it is evident that this effective background 
has an apparent future event horizon only for n > 1 and 
we can trust the decoupling of long wavelength modes 
from short wavelength modes. On the other hand, for 
n < —2, we can use the symmetry of S"/S = n(n+ 1)/t 2 
under n — > — (n+ 1) to obtain a consistency relation for 
those values, once we have a consistency relation valid 
for n > 1. 

It is interesting to compare our result in (|20[) . with 
the calculations of [ill and [13j. In [13J, the comoving 
curvature perturbation was replaced with the perturba- 
tion of a test scalar field. On the other hand, if we use 
the relation £ = S(j>/y/2e, the dominant interaction of 
the metric perturbations with the electromagnetic field 
can be obtained by Taylor expanding the coupling in 
the infiaton fluctuations as A(</)) = A(</> c ) + d^X(4> c )5(j). 
In this way, one obtains the linear coupling between 
the inflation fluctuation and the electromagnetic field 
£-5<pBB = — (l/4:)d l fMd))S(f)F 2 . This is same as the cou- 
pling considered in }13j | and we are able to reproduce their 
results in the squeezed limit when using AC = —d^XHScj) 
in order to write the consistency relation as 

(tyfokOBfa.k-O.BCTr.ka)) 

= ^( 27r ) 3 ( 5( 3 )(k 1 + k 2 + k^Ps^PBih) (21) 

and inserting the specific form of the coupling X(4>) = 
exp(2(f>/M) used there. The power spectrum of the 
scalar field fluctuations is defined as (S(/)(t, k)(50(r, k')) = 
(2t:) 3 S^ (k + k')Ps < p(k). By means of numerical checks, 
it was observed that in the squeezed limit, the final re- 
sults are independent of n (l3l |. We can now understand 
straightforwardly from (|2"TT) why this has to be the case. 
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In [11J, the correlation function in (1201) was calculated 
for the actual metric perturbation. The form of the ac- 
tion that they used is slightly more complicated than 
C = — (1/4)<90A(</>)<50F 2 to the leading order in slow- 
roll. However, it can be shown that the extra terms 
in their interaction Hamiltonian are proportional to a 
total derivative 21| and in fact, using a full gauge fix- 
ing in the Arnowitt-Deser-Misner formalism [22|, it was 
shown in Q that at linear order, the inflaton fluctuation 
5<f> is related to the comoving curvature perturbation by 
C = S(f)/^/2e and the relevant part of the action to lead- 
ing order in slow-roll is indeed C = —(l/£)d<p\((j))5(l)F 2 . 
Therefore, the squeezed limit result of should agree 
with that of 
what 



Foundation. 
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13j, when using £ = S(f>/^/2e. This is not 



finds which indicates that there is numerical 
factor wrong in the results of llll . We have therefore 
carefully checked the results of [11| and a detailed calcu- 
lation shows that the correct calculation agrees with the 
results of [HI and with the consistency relation in the 



squeezed limit [21 



We have thus shown in (|20|) that in a generic model 
where magnetic fields are generated during inflation, 
there exists a new type of consistency relation for the 
three point cross correlation function of the comoving 
curvature perturbation with two powers of the magnetic 
field in the squeezed limit where the momentum of the 
curvature perturbation vanishes. Moreover, this consis- 
tency relation can also be written for the correlation of in- 
flaton field fluctuations with the magnetic field as in ([2~Tj) . 
Similarly, when the electromagnetic field is coupled to a 
curvaton-like scalar field (instead of the inflaton) ignor- 
ing gravity in a fixed de Sitter background, the correla- 
tion of the field fluctuation with the magnetic field in the 
squeezed limit is again given by (|21l) . 

Although the underlying model suffers from the strong- 
coupling problem 0, Qi| , it is an archetypical model of 
primordial magnetogenesis which has been widely stud- 
ied in the literature. It will be interesting to see if the 
strong coupling problem can be avoided (perhaps along 
the lines suggested in [HI), an d also, if the results ob- 
tained here can be applied to other models. For instance, 
the gauge field might instead play the role of a vector cur- 
vaton leading to some amount of statistical anisotropy 
and anisotropic non-Gaussianity though it requires con- 
siderable fine tuning of the gauge coupling [23 1. 
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